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ABSTRACT 

Searches for slow radio transients and variables have generally focused on extragalactic populations, 
and the basic parameters of Galactic populations remain poorly characterized. We present a large 
3 GHz survey performed with the Allen Telescope Array (ATA) that aims to improve this situation: 
ASGARD, the ATA Survey of Galactic Radio Dynamism. ASGARD observations spanned 2 years 
with weekly visits to 23 deg^ in two fields in the Galactic Plane, totaling 900 hr of integration time on 
science fields and making it significantly larger than previous efforts. The typical blind unresolved 
source detection limit was 10 mJy. We describe the observations and data analysis techniques in 
detail, demonstrating our ability to create accurate wide-field images while effectively modeling and 
subtracting large-scale radio emission, allowing standard transient-and- variability analysis techniques to 
be used. We present early results from the analysis of two pointings: one centered on the microquasar 
Cygnus X-3 and one overlapping the Kepler field of view {£ = 76°, b = +13.5°). Our results include 
images, catalog statistics, completeness functions, variability measurements, and a transient search. 
Out of 134 sources detected in these pointings, the only compellingly variable one is Cygnus X-3, 
and no transients are detected. We estimate number counts for potential Galactic radio transients 
and compare our current limits to previous work and our projection for the fully-analyzed ASGARD 
dataset. 

Subject headings: radio continuum: general — surveys — techniques: interferometric 



1. INTRODUCTION 

The technological developments of the past few decades 
have led to an explosion of interest in the astronomical 
time domain. There has been a recent boom at radio 
wavelengths, where the survey capabilities of new and 
upgraded facilities represent order-of-magnitude improve- 
ments over their predecessors (e.g., Perley et al. 2011; 
Verheijen et al. 2008; Giovanelli et al. 2005; DeBoer et al. 
2009; Jonas 2009; Carilh & Rawlings 2004). Among 
the motivators for the construction of these facilities are 
known or strongly-supported classes of highly-variable 
slow extragalactic radio emitters such as active galactic 
nuclei (Lister et al. 2011), orphan 7-ray burst afterglows 
(Frail et al. 2001; Levinson et al. 2002), radio supernovae 
(Weiler et al. 2002; Brunthaler et al. 2009; Muxlow et al. 
2010), and tidal disruption events (Rees 1988; Bloom 
et al. 2011; Berger et al. 2012). We follow other authors 
in defining slow variables as those whose emission evolves 
on timescales ^1 s; this approximately corresponds to 
those that emit via incoherent, rather than coherent, pro- 
cesses, and are typically identified using image-domain 
techniques. The relatively sparse prior exploration of the 
dynamic radio sky additionally highlights it as an exciting 
discovery space (Cordes et al. 2004; Fender & Bell 2011; 
Frail et al. 2012). 

Much work has recently gone into the characterization 
of the population of highly- variable slow extragalactic 
sources. Archival studies have used existing large- area 
surveys, usually observed at 1.4 or 4.9 GHz, to search 
for rare events (Levinson et al. 2002; Bower et al. 2007; 
Beh et al. 2011; Bower & Saul 2011; Ofek & Frail 2011; 
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Bannister et al. 201 la, b; Thyagarajan et al. 2011, 2012; 
Frail et al. 2012). Followup observations of candidate 
or confirmed highly-variable sources have been used to 
characterize the detailed properties of individual objects 
(de Vries et al. 2004; Gal- Yam et al. 2006; Muxlow et al. 
2010; Berger et al. 2012). Finally, dedicated surveys 
have been devoted to the systematic discovery of highly 
variable radio sources. Several have been undertaken with 
the Allen Telescope Array (ATA; Welch et al. 2009), one 
of the first radio observatories exphcitly designed to be an 
efficient survey instrument, including the ATA Twenty- 
centimeter Survey (Croft et al. 2010a,b, 2011) and the 
Pi GHz Sky Survey (Bower et al. 2010, 2011, PiGSS-I 
and PiGSS-II hereafter; Croft et al, 2012, in prep.). Ofek 
et al. (2011) describe most of these previous studies in 
greater depth. Frail et al. (2012) summarize the state of 
the field and find that reliable detection of these sources 
remains both a challenge and an opportunity. 

In contrast, there has been relatively little work to 
investigate the population of highly- variable slow Galac- 
tic radio emitters. Source classes contributing to this 
population include X-ray binaries (Waltman et al. 1995; 
Marscher & Brown 1975), active stellar binaries (Hall 
1976; Bopp & Fekel 1977; Eker et al. 2008), cool dwarfs 
(Berger et al. 2001; Berger 2002; Hallinan et al. 2007), 
pre-main-sequence stars (Bower et al. 2003; Forbrich 
et al. 2008; Salter et al. 2008), and fiare stars (Giidel 
2002; Jackson et al. 1989; Richards et al. 2003). To- 
wards the Galactic Center (GC), several intriguing sources 
of ambiguous nature have been discovered, including 
the Galactic Center Transient at 1.36 GHz (Zhao et al. 
1992), A1742-28 at 0.96 GHz (Davies et al. 1976), 
CXOGC J174540.0-290031 at a variety of frequencies 
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(Bower et al. 2005), GCRT J1745-3009 at 0.33 GHz 
(Hyman et al. 2005, 2006), and GCRT J1742-3001 at 
0.235 GHz (Hyman et al. 2009). An overall increase in 
the prevalence of apparent radio variability is expected to- 
wards the Galactic plane (GP), and the GC in particular, 
due to interstellar scintillation (e.g., Spangler et al. 1989; 
Rickett 1990; Ghosh & Rau 1992; Gaensler & Hunstead 
2000; Loveh et al. 2008; Ofek & Frail 2011). Pulsars 
and several other well-known classes of variable Galactic 
sources are not discussed here because they do not fall 
into the "slow" category; see Cordes et al. (2004) and 
references therein. 

For many years, the best available data on Galactic 
radio variability came from the work of Taylor & Gregory 
(1983) and Gregory & Taylor (1986), who used the NRAO 
91-m transit telescope to repeatedly survey the GP at 
5 GHz. Recent work has begun to expand and update 
these results. Hyman et al. (2002, 2003, 2006) have peri- 
odically monitored the GC at long wavelengths (0.3 GHz) 
with the VLA and GMRT, and as indicated above they 
have discovered several unusual transient sources. Becker 
et al. (2010) conducted an archival search for GP tran- 
sients at 5 GHz by comparing the VLA survey described 
in Becker et al. (1994) with CORNISH, the "Co-Ordinated 
Radio 'N' Infrared Survey for High-mass star formation" 
(Purcell et al. 2008), discovering a population of highly 
variable sources and analyzing their statistical properties. 
Most recently, Ofek et al. (2011) conducted a 5 GHz VLA 
variability survey at low Galactic latitudes with rapid 
multiwavelength followup, discovering a candidate tran- 
sient source and measuring more statistical properties of 
the population. 

In this paper we present ASGARD: the ATA Survey 
of Galactic Radio Dynamism. Its primary goals are to 
perform a deep search for Galactic radio transients and 
to measure the variability properties of a wide variety 
of Galactic radio sources on day-to-year timescales. To 
this end, we repeatedly observed 24 pointings near the 
GP at 3 GHz with the ATA, visiting most pointings on a 
^weekly cadence over the course of ~2 yr and obtaining 
~900 hr of integration on our science pointings. The 
large field of view (FOV) of the ATA allows us to cover 
a relatively large region on the sky, and our frequent 
visits provide thorough sampling of variability on a range 
of timescales. The compact configuration of the ATA 
also makes our observations sensitive to static extended 
structures in the GP, such as nonthermal radio filaments 
(Yusef-Zadeh et al. 1984; Law et al. 2008), H II regions 
(Brogan et al. 2003; Nord et al. 2006), and supernova 
remnants (Gray 1994). 

We proceed by describing the ASGARD observations 
(§2) and data processing (§3) in detail. This work con- 
cerns itself with a subset of the whole dataset (§3.1), 
representing ~10% of the expected usable observations, 
while another portion of the dataset has already been 
described elsewhere (Williams et al. 2011). We present 
first results (§4) derived from analysis of this subset, in- 
cluding deep images, source catalogs, a transient search, 
and variability statistics. Finally we discuss our current 
conclusions and the prospects of the fully-analyzed survey 
(§5), populating a "log- A/"/ log- 5*" plot for Galactic radio 
transients (Figure 12). 

2. OBSERVATIONS 



We used the ATA to monitor two fields: an area around 
the Galactic center (GC) spanning -4.5° < f < 8°, |6| < 
2°; and ~5 deg^ towards Cygnus including the highly 
radio- variable microquasar Cyg X-3 (Gregory & Kronberg 
1972) and a portion of the FOV of the Kepler (Koch et al. 
2010) mission. Coordinates of the ASGARD pointing 
centers are listed in Table 1, along with the observing 
time devoted to each pointing over the course of the 
survey. Throughout this work we use the word "field" to 
refer to either of the two survey regions (GC and Cygnus) 
and the word "pointing" to refer to a specific pointing 
center that was observed. 

While some ASGARD observations were conducted 
with complete control of the ATA, the vast majority 
were conducted commensally with a SETI (Search for 
Extraterrestrial Intelligence) survey (Blair & ATA Team 
2009). Both of these projects were designed to increase 
the likelihood of discovering rare Galactic events by tar- 
geting regions with high source densities, as well as to 
take advantage of the data stream provided by the Kepler 
mission, which motivated the basic choice of survey fields. 
The division of the survey area into two fields was also a 
practical choice given year-round observations, since the 
GC region is seasonally difficult or impossible to observe 
from the ATA site (latitude +40.8°). The Cygnus region, 
on the other hand, can be observed year-round, also has 
a very high source density, and contains the benchmark 
source Cyg X-3. The observed GC field is as centered on 
Sgr A* as possible given visibility constraints. The overall 
footprint of each field was determined by the joint needs 
of the commensal SETI search and ASGARD. The former 
required a certain dwell time per pointing to sequentially 
survey targets using the ATA digital beamformer back- 
ends. The latter aimed for a weekly revisitation cadence. 
Combining these requirements with the typical weekly 
survey time allocation determined the total footprint that 
could be observed. The system used to organize and exe- 
cute the commensal observations is described in Williams 
(2012). 

The pointings in the GC field fall on an 11° x 2° square 
grid in Galactic coordinates with a spacing of 1° and a 
northeast corner located at ^ = 6.5°, h = 0.5°. The half- 
power beam width (HPBW) of the ATA is approximately 
(3.5 deg GHz)/i/obs (Huh et al. 2010; Harp et al. 2011, 
but see §3.5), so that the pointing centers range from 
slightly oversampled to critically sampled for ASGARD 
observing frequencies of ~l-3 GHz. The Cygnus field 
consists of an analogous 2x2 grid, with Cyg X-3 located 
at the northwest pointing center {£ = 79.8°, b = 0.7°), as 
well as a disjoint pointing toward £ = 75.8°, b = 13.5° 
{a = 19^21^^24% S = 44°00'00'^ ICRS J2000), which lies 
within the Kepler mission FOV. The total footprint of 
the survey on the sky is ~23 deg^ if each pointing is taken 
to image a circle with a diameter of the nominal HPBW. 

At the time of the observations the ATA frontend con- 
sisted of forty- two 6.1-m offset Gregorian dishes. The 
backends used for this work were two FX correlators with 
bandwidths of 104.9 MHz divided into 1024 channels and 
a dump time of 10 s. The correlators accepted 64 "antpol" 
inputs, meaning that they could perform full- Stokes cor- 
relation of 32 antennas or, hypothetically, single-Stokes 
correlation of 64 antennas. Because of the desirability 
of full- Stokes coverage and an ongoing program of feed 
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retrofits, the correlators generally accepted data from ~32 
distinct antennas. The set of antennas used varied with 
time due to maintenance or hardware failures, and the 
two correlators did not necessarily have identical antpol 
inputs. 

Each observing session ("epoch") generally began with a 
long (~30 min) observation of a bright, unresolved source, 
usually one of 3C48, 3C 147, 3C 286, or 3C295. These 
observations were necessary for delay calibration of the 
beamformers used in the SETI survey but also provided 
excellent bandpass and flux density scale calibration data 
for ASGARD. Science pointings were observed with peri- 
odic (every ~45 min) visits to a nearby phase calibrator. 
Because of the low dechnations of the GC pointings, these 
were generally observed at relatively low elevations with 
constrained hour angle ranges. 

The ATA feeds have a high-bandwidth design using 
a log-periodic architecture, where the location of the 
active region on the feed varies with the observing fre- 
quency. Because the optics of the ATA reflectors are 
frequency-independent, each feed is mounted on a piston 
drive so that the appropriate part of the feed may be 
moved to the optical focus for each observation. Focus 
positions are identifled by the optimal corresponding ob- 
serving frequency, i.e., a focus position of i^foc = 1.4 GHz 
provides the best conflguration for observations at that 
frequency. The relationship between the piston position 
and observing frequency is based on both theoretical and 
empirical analysis (Harp et al. 2011). Defocused observa- 
tions are possible with some loss in system performance. 
For z/foc ^ O.9i^obs5 the penalty is slight, while for lower 
focus frequencies (active region too close to the secondary) 
sensitivity degrades and the primary beam (PB) broadens. 
For ASGARD observations made in complete control of 
the array, the focus position was set optimally, usually at 
3.14 GHz. For commensal observations, the focus position 
was set at 1.90 GHz. 

The overall coverage statistics of ASGARD observations 
are recorded in Table 2. We divide the observations 
into four campaigns: three seasons of GC observations 
during the summers of 2009-2011, and regular Cygnus 
observations from 2009 November to 2011 April. For 
reasons discussed in the following section, we isolate the 
set of observations made with the ATA correlators tuned 
to sky frequencies of 3.04 and 3.14 GHz, which we refer 
to as the "3 GHz" observations. These are the same 
frequencies used by the PiGSS survey. The overall amount 
of observatory time dedicated to the project was 1650 hr, 
with 902 hr spent observing science targets. 

3. ANALYSIS 

We mostly use standard techniques to calibrate and im- 
age the data, using tools from the MIRIAD data reduction 
package (Sault et al. 1995) for many steps, with several 
customized steps implemented in the Python^ program- 
ming language via the package miriad- python (Williams 
et al. 2012). We use our images to construct a catalog 
of every compact source detected in the survey and to 
measure flux densities (or upper limits) for every observa- 
tion of every source. We perform a variability analysis on 
this photometric dataset. We deflne a criterion by which 
sources may be classified as "transient" or not, but the 

^ http://python.org/ 



classification does not affect the variability analysis, and 
no transient sources are detected. 

The chief difficulty in the particular case of ASGARD is 
successful imaging of the significant large-scale structure 
(LSS) present in virtually all ASGARD fields of view, with 
the notable exception of the Kepler pointing. Not only 
is LSS imaging generally challenging, but by the nature 
of the survey most ASGARD epochs have sparse hour 
angle coverage. Fortunately, the LSS emission is expected 
to be time-invariant, while any astrophysical transients 
will be unresolved by the ATA (resolution at 3 GHz). 
All observations of a given pointing can therefore be 
combined to produce an LSS model that can then be 
subtracted in the visibility domain from each epoch's 
observations. To ease the measurement of the variability 
of compact sources, we exclude these sources from the 
LSS model, causing them to remain in the per-epoch 
images. Putting aside for the moment the important 
issue of subtraction errors, this approach yields per-epoch 
images consisting only of compact sources upon which 
standard transient-search techniques may be used. We 
describe our implementation of this general approach 
below. 

3.1. Subset of Data Presented in This Work 

The first season of GC observations was performed at 
sky frequencies of 1.43 and 2.01 GHz. It was found, how- 
ever, that there was substantial broadband interference at 
these frequencies that presented numerous challenges for 
data analysis (Williams 2010). Subsequent observations 
primarily used the PiGSS 3 GHz setup, although a few 
observations were made at other frequencies. The analysis 
presented in this work is restricted to 3 GHz observations. 
In Tables 1 and 2 we provide coverage statistics as com- 
puted for the 3 GHz subset of the complete ASGARD 
dataset. 

In this work, we analyze two particular ASGARD point- 
ings. The first is towards the highly-variable source 
Cyg X-3, which allows us to demonstrate the detection of 
variable radio sources embedded in complex, large-scale 
emission. The second is the Kepler pointing, which allows 
us to investigate the performance of our techniques in 
a field without unusual imaging challenges. We demon- 
strate the analysis of 29 epochs of the Cyg X-3 pointing, 
representing 34% of the 3 GHz Cyg X-3 dataset by num- 
ber of epochs and 17% by raw data volume, and six epochs 
of the Kepler pointing, five of which also involved visits 
to the Cyg X-3 pointing. These epochs nearly completely 
sample the timespan of the 3 GHz Cygnus campaign, 
ranging from 2010 February 03 to 2011 April 04, and are 
roughly uniformly spaced in time. We have processed 
and imaged epochs surrounding the Cyg X-3 radio fiares 
of 2010 May (Bulgarelh et al. 2010) and 2011 March 
(Kotani et al. 2011), but do not run our full pipeline on 
the latter epochs, of which there are six, because the 20 Jy 
fiare leads to severe imaging problems related to dynamic 
range limitations. Tests of our transient detection process 
confirm, however, that it succeeds in this trivial case. 
Other work describes the ATA-42 observations of these 
events in more detail (Williams et al. 2011). A somewhat 
larger portion of our observations has been processed and 
analyzed but without significant multi-epoch coverage of 
other pointings, and so is not presented in this work to 
maintain a clear focus on the better-covered Cyg X-3 and 
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Kepler pointings. 

3.2. Calibration & Flagging 

Radiofrequency interference (RFI) was an intermittent 
problem during the 3 GHz observations. The forms of RFI 
most commonly encountered were narrow-bandwidth (1-5 
channel) tones that affected the majority of baselines with 
a ^100% duty cycle, wider-bandwidth tones (~5 MHz) 
that affected moderate numbers of baselines with a mod- 
erate duty cycle, and brief (< 1 dump) broadband bursts 
affecting all baselines. Problems in the digital hardware 
(e.g., overheating) could also cause RFI-like effects, typi- 
cally manifesting themselves as complete corruption of one 
half or one quarter of the spectrum for certain baselines. 
In the results we report, RFI was primarily excised from 
the data manually, using standard MIRIAD tools and 
an interactive, graphical visibility visualizer for the RFI 
with more complex time/frequency structure (Williams 
et al. 2012). Several approaches to automatic flagging 
have also been pursued with ATA data (e.g., Keating 
et al. 2010; Bower et al. 2010) and their integration into 
the ASGARD pipeline is being investigated. 

Standard bandpass and gain calibration techniques are 
used. The long calibration observations at the beginning 
of each epoch are used to set the flux density scale, refer- 
encing to Baars et al. (1977). For the GC observations, the 
gain calibrator was NRAO 530 {I = 12.03, 6 = +10.81), 
while for the Cyg X-3 observations it was usually BL Lac. 
Because this latter source is variable and usually ~10% 
linearly polarized, gain parameters for it are derived from 
the bandpass observations whenever possible, treating X 
and Y feeds separately. Our observing program did not 
allow for planned observations of polarimetric calibrators 
over wide hour angle ranges, so we are unable to solve 
for the frequency-dependent leakage terms that would 
be required for polarimetric calibration of the ATA (Law 
et al. 2011b) on an epoch- by-epoch basis. 

3.3. Imaging and Source Extraction 

After calibration, the data are averaged down to 16 
spectral channels of 6.5536 MHz bandwidth per correlator 
and are converted to CASA (McMulhn et al. 2007) format 
for imaging. This is necessary because MIRIAD does not 
implement any wide- field imaging algorithms, which we 
have found to be necessary for our analysis. In particular, 
without the use of techniques such as polyhedral imaging 
(Cornwell & Perley 1992) or iL'-projection (Cornwell et al. 
2008), sources far away (^0.8°) from phase center do not 
deconvolve well and acquire an hour angle dependence 
in their position, which is severely problematic for both 
LSS modeling and photometric extraction. In the steps 
we describe below, the averaged channels are gridded 
using multifrequency synthesis (Sault & Wieringa 1994) 
and imaged using K;-projection with 128 planes (but not 
polyhedral techniques). Images are 2048x2048 with a 
pixel size of 10^' and thus span approximately five times 
the HPBW at 3 GHz. 

We first construct a deep image for each pointing. Be- 
cause these images are used to generate several important 
data products, the imaging techniques vary from pointing 
to pointing depending on what produces the best results 
empirically. For the Kepler pointing, we use the Cotton- 
Schwab deconvolution method (Schwab 1984; Cornwell 



et al. 1999), while for Cyg X-3 we use the maximum- 
entropy algorithm (Gull & Daniell 1978) with a Gaussian 
prior. To achieve consistent LSS sampling without requir- 
ing detailed primary beam models, the data contributing 
to each deep image have a consistent ATA feed focus 
position. Imaging artifacts, rather than thermal noise, 
currently limit the image quality, so the use of only part 
of our data to form the deep images does not significantly 
affect their sensitivities. 

The deep images are used to construct a catalog of com- 
pact sources, the properties of which are described in §4.2. 
Sources are detected using a combination of the MIRIAD 
task sfind and manual inspection to check for missed de- 
tections and reject dubious ones. (Given the resolution, 
sensitivity, and footprint of ASGARD, these techniques 
are scalable to the whole survey.) Source positions, shapes, 
and mean flux densities are cataloged. Because the pri- 
mary aim of ASGARD is to study variability, our catalog 
does not initially include sources that are marginally de- 
tected in the deep images, since such sources wih typically 
not be detectable in the epoch images. The epoch im- 
ages, however, are searched for uncataloged sources as 
described below, so that any variable source that becomes 
detectable during an epoch will eventually be included in 
the catalog. 

Each ASGARD pointing is also associated with a static 
LSS model. In the Kepler pointing, this model is blank; 
for the rest, the model is derived from the deep image and 
the compact source catalog, using either deconvolution 
of a variant deep image in which the compact sources 
have been subtracted, or source- fitting techniques on the 
model image. The latter approach can be helpful be- 
cause maximum-entropy deconvolution tends to model 
unresolved sources as Gaussians about the size of the 
synthesized beam. Because there are substantial numbers 
of Cyg X-3 observations made at focus positions of both 
1.90 and 3.14 GHz, we generate one LSS model of this 
pointing for each focus position, so that the approxima- 
tions of our primary beam modeling scheme (§3.5) do not 
lead to avoidable LSS subtraction errors. 

Images from individual epochs are made the same way 
as the deep images, except that before imaging the appro- 
priate LSS model is subtracted from the u-v data, and 
baselines shorter than 50 m are not imaged. Deconvo- 
lution of the individual epoch images is performed with 
800 iterations of CASA's "wide-field" implementation of 
the Hogbom CLEAN algorithm (Hogbom 1974). The 
restoring beam is chosen automatically by the imaging 
software and has a typical size of 80^' x AO" . In two epochs 
of Cyg X-3 observations, there are small but discernable 
errors in the flux density scale that lead to noticeable 
residuals in the LSS subtraction. We fixed the scales in 
these epochs by trial- and-error, adjusting and reimaging 
the data and visually assessing the magnitude of the LSS 
subtraction residuals. The correction factors are -1% and 
5%, and we estimate that they are uncertain by about 1%. 
(Because the LSS is strong and unvarying, it is possible 
to be relatively precise.) We discuss our investigations 
into rigorous, global cross-calibration of the flux density 
scale in §3.6. 

Lightcurves of the cataloged sources are derived using 
image-domain fitting on the LSS-subtracted individual 
epoch images. Each fit holds the source position and 
shape fixed but allows the total flux density to vary. Very 
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close sources are fit simultaneously. Because the reality, 
positions and shapes of these sources are known from the 
deep image, we use a fairly weak constraint and consider 
a source to be detected if its fitted flux is more than three 
times the local background rms. Undetected sources are 
cataloged with an upper limit of this theshold. 

During the source fitting process, a residual image is 
generated by subtracting off the best-fit image-domain 
model of every detected source. There are generally sub- 
traction residuals around each source of peak magnitude 
<10% of the source flux, with both positive and negative 
components. The mean residual is another ~10% smaller 
because the fitting process tends to minimize this value. 
Allowing the source shapes and/or positions to vary yields 
smaller residuals by definition, but the derived parameter 
values vary much more from epoch to epoch than instru- 
mental and observational uncertainties lead us to expect. 
We instead attribute the subtraction residuals primarily 
to u-v calibration errors. In this case our choice to fix 
the source positions and shapes gives a more realistic 
assessment of flux uncertainties by avoiding overfitting of 
the data. 

3.4. Detection of Uncataloged Sources 

We use sfind to search for any residual sources in the 
individual epoch images after subtraction of both the 
LSS and the cataloged compact sources. The LSS is 
subtracted in the u-v domain but the compact sources 
are subtracted in the image domain. In our analysis, 
a "transient" is any source detectable in an individual 
epoch image that is not detectable in the corresponding 
deep image. Sources discovered via sfind are added to 
the ASGARD catalog and so are subsequently processed 
in the same way as all others. For those pointings in 
which the LSS model is derived from a deep image, such 
a source's contribution to the deep image propagates into 
the LSS model. This mean flux density is equal to the 
bright-epoch flux density diluted by ~l/n, where n is the 
number of epochs. This contribution is subtracted from 
the per-epoch images used for residual source detection 
but will not significantly alter detectability for n more 
than a few. 

Previous searches for radio transients have typically 
been dominated by false positives (Frail et al. 2012). Our 
procedure involves multiple rounds of sky modeling and 
subtraction which will inevitably leave artifacts as well. 
Our transient-detection step therefore has stringent detec- 
tion limits and cross-checks for systematic effects. We use 
the false-discovery rate (FDR) algorithm in sfind (Hop- 
kins et al. 2002). The background rms is computed in 
64 X 64-pixel boxes (sfind keyword rmsbox) and the target 
FDR is set to 0.5% (sfind keyword alpha). To assist in 
computing detection limits, sfind was modified to report 
an estimated minimum detectable source flux density in 
the event that no sources were found in an image, basing 
this value on an estimated FDR "p- value" that would be 
needed to have yielded a source detection. The source 
shapes reported by sfind are deconvolved from the syn- 
thesized beam. (Recall that a typical synthesized beam 
size in our images is 80^' x 40" .) Sources smaller than 
the synthesized beam have their shape fixed to that of 
the beam and have their parameters refit (sfind option 
psfsize). Those that are otherwise incompatible with the 
synthesized beam shape (e.g., both elongated with per- 



pendicular major axes) are also treated as point sources. 
This is a conservative approach because, as described in 
the next paragraph, we reject extended sources. 

Sources reported by sfind are filtered according to sev- 
eral criteria. Sources for which sfind does not report a 
positional uncertainty, usually indicative of a very poor 
fit, are rejected. Because genuine astrophysical transients 
will not be resolved by the ATA, sources in which the 
product of the deconvolved major and minor axes exceeds 
7000 arcsec^ are rejected, as are those in which the decon- 
volved major axis exceeds 130 arcsec. Sources for which 
the modeled primary beam attenuation exceeds 98.9% 
(separation of 3crpB, where apB = HPBW/V81og2) are 
rejected, although our later analysis uses much more con- 
servative PB attenuation cutoffs. Finally, newly-detected 
sources near previously cataloged steady sources are also 
rejected. The match radius for this test is 50/ arcsec 
where / depends on the cataloged total flux density Sk 
of the known source. For Sk < 20 mJy, / = 1; for 
Sk > 2.97 Jy, / = 10; and for intermediate values, 
/ = 21og(5'/c/20 mJy). All of the above cutoff values 
were determined by examining the properties of the tran- 
sient candidates that were both significantly detected and 
obviously spurious. Remaining candidates are examined 
manually as described in §4.4. 

3.5. Primary Beam Modeling 

In order to compute accurate flux densities, sensitivities, 
and sky models, we must account for the primary beam 
of the ATA. Hull et al. (2010) analyzed the primary 
beam of the ATA using data from the PiGSS survey, 
finding HPBW = 1.10 ib 0.01° for a circular Gaussian 
primary beam model. Because many of our observations 
are performed with the same frequency configuration as 
PiGSS, we adopt this value assuming a frequency 
dependence, i.e. HPBWoptimai = (3.40 deg GHz)/i^obs 
for a mean PiGSS observing frequency of 3.09 GHz. The 
numerator is slightly smaller than the more generic value 
of 3.5 deg GHz reported by Harp et al. (2011). 

Our PB modeling is complicated, however, by the fact 
that a substantial number of 3 GHz observations are made 
with the focus set to 1.90 GHz. Although we do not have 
observations specifically aimed at characterizing the PB 
shape of 3 GHz observations at this focus setting, we 
measure this value from the data in two ways. Firstly, 
our observations of the Cyg X-3 pointing have extensive 
coverage in both this focus position and in the optimally 
focused position. By comparing compact source flux 
densities in deep images made from the two sets of data, 
we obtain HPBWigoo = 3.92 deg GHz/i^obs- We also 
compare the ATA-apparent flux densities of point sources 
observed in the Kepler pointing to those from the NVSS 
catalog (see §4.2), assuming a typical source spectral index 
a of -0.7 {Si, oc i^"), and find a factor of 3.96 deg GHz. We 
use the former value. Below, we denote the PB correction 
factor as applied to a particular location /pb, defined 
such that /pB > 1. 

Holographic measurements of the ATA dishes suggest 
that the PB pattern becomes increasingly noncircular as 
the focus moves away from its optimal setting (Harp et al. 
2011). Due to alt-az mount of the ATA dishes, the PB ro- 
tates on the sky over the course of an observation. Proper 
accounting for such an effect would need to occur during 
the Fourier inversion process, which neither MIRIAD nor 
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the standard CASA imager are capable of doing. We are 
thus unable to measure, or compensate for, this effect. 
(PB rotation can be dealt with approximately by imaging 
the data in blocks of similar hour angle, but this reduces 
u-v coverage and thus hampers the deconvolution of our 
complex fields.) The measurements of Harp et al. (2011) 
suggest that in our configuration the axial ratio should 
be about 10%. 

We also assume that each dish has an identical PB 
pattern and that all dishes are pointed identically. These 
assumptions allow us to handle PB correction in the image 
domain and are relied on in some of the analysis that 
follows. The ability to relax these assumptions, like the 
ability to model PB rotation, relies on either painstaking 
iterative approaches (Wright & Corder 2008) or support in 
imaging software that is not yet widespread; two notable 
implementations are the MeqTrees system (Noordam & 
Smirnov 2010) and a derivative of the CASA imager 
equipped with the A-projection algorithm (Bhatnagar 
et al. 2008). Adapting our analysis to a pipehne in which 
PB correction occurred during the imaging process would 
require actual measurement, rather than calculation, of 
the spatial variation in image noise, but would not require 
major conceptual changes. 

3.6. Multi-Epoch Photometry 

Our imaging and source extraction pipehne yields multi- 
epoch photometry for our catalog sources, though for 
some faint sources our measurements yield mostly upper 
limits. Each source additionally has a flux density mea- 
surement from the deep image in which it is best detected. 
Some bright sources are detectable in the deep images of 
multiple pointings. 

Our images are not corrected for the attenuation of the 
ATA primary beam, so we must correct our flux density 
measurements for this effect. In the subset of data we 
consider, all Kepler observations are made with a focus 
setting of 1.90 GHz, and so the primary beam correction 
for a given source is the same in every epoch. The Cyg X-3 
observations are made with focus settings of both 1.90 
and 3.14 GHz, so the primary beam correction can vary 
from epoch to epoch. The deep Cyg X-3 image derives 
from the optimally-focused data. Comparison of the 
Cyg X-3 flux densities is complicated by the fact that we 
use different LSS models for each focus setting. While this 
is necessary to accurately subtract LSS from each epoch 
as well as possible, differences in the LSS model around 
each sourch can lead to an additive flux density offset 
beween measurements made at different focus settings, 
above and beyond any multiplicative errors caused by the 
limitations of our analytic primary beam models. In this 
work, we choose to consider only measurements from a 
consistent focus setting, choosing the one that resulted in 
the higher mean detection significance. For sources close 
to the pointing center, these are generally the optimally- 
focused data, whereas for sources far from the pointing 
center, the broad defocused primary beam results in more 
significant detections. 

Some images yield particularly poor photometric results 
and we discard these measurements. These are generally 
either epochs with very few contributing data or those 
from 2011 March in which Cyg X-3 was undergoing a 
major (~20 Jy) fiare, leading to significant dynamic range 
issues. Even without the inclusion of the latter flare in 



our analysis, Cyg X-3 is the most signiflcantly variable 
source in our dataset. 

We set the uncertainty on each flux density measure- 
ment to be 

= ^rms + (0.055)2 ^ 5 ^jy)2^ 

where arms is the background rms of the relevant image- 
domain fit and S is the measured value. The constants in 
this equation, which are in the range typically encountered 
in radio interferometry, were chosen to yield a plausible 
distribution of Pc values (see §4.5) for the least-variable 
survey sources. We do not include the uncertainty of the 
fiux density scale correction factor for the two epochs that 
were rescaled (§3.3) or attempt to quantify the uncertainty 
on the PB correction. We investigated a correction for 
CLEAN bias (Becker et al. 1995; Condon et al. 1998) 
but did not find compelling evidence that this improved 
our results, so we do not include such a correction in our 
analysis. We also investigated but did not use parabolic 
photometric models as employed for some sources by 
Bannister et al. (2011a,b). The full ASGARD dataset 
should have enough flux density measurements to be able 
to investigate the impact of these and other techniques 
in a more statistically rigorous manner. 

Bannister et al. (2011a,b) and Ofek et al. (2011) used 
"post-imaging calibration" (PIC) techniques in which they 
cross-calibrated their photometry assuming that most 
sources do not vary systematically and that each image 
is subject to a multiphcative flux density scale correction. 
The flux density scale correction that we describe in §3.3 
contrasts in that it happens before the imaging process, 
is only applied to two problematic epochs, and uses cali- 
bration factors determined in a clearly more ad-hoc way. 
We investigated but did not use a "pre-imaging calibra- 
tion," in which we correct the flux density scale of every 
epoch's visibilities using a correction factors derived from 
the photometric data as in PIC. Besides being extremely 
computationally costly, since applying the calibration 
requires reimaging the entire survey, we found in small- 
scale tests that the calibration factors were not stable 
over multiple iterations of the calibration. We also found 
that this approach did not correct the LSS subtraction 
residuals in the two problematic epochs nearly as well as 
we could manually. We additionally considered but did 
not use a true PIC, operating only in the photometric 
domain. As was also found in PiGSS-II, we found that 
such a calibration sometimes had encouraging results, 
especially for bright sources, but that its effect on the 
photometry of faint sources appeared to range between 
neutral and negative. As such we decided to err on the 
side of simplicity and fewer data transformations. 

For those compact sources that are found within a 
region of more extended emission, our LSS subtraction 
technique naturally leads to the possibility of an additive 
flux density bias in our measurements. We constructed a 
set of images in which the LSS model was added back to 
each epoch image after convolution with the synthesized 
beam. We then extracted photometry again using a 
specialized routine that simultaneously flt for the flux 
density of a source in each image on top of a local, constant 
background term. This more elaborate procedure did not 
produce more consistent results, and we found that our 
LSS subtraction procedure did not lead to consistent flux 
density biases. 
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In Figure 1 we compare our deep and epoch flux density 
measurements after correcting the flux density scale of 
the two epochs that needed it, accounting for PB atten- 
uation, discarding low-quality points, and augmenting 
the uncertainties as per Equation 1. There is very good 
agreement in the scales between the two. Cyg X-3 is 
clearly the most significantly variable source we observe. 
The nominally very bright sources are all detected at very 
large distances from the pointing center and are subject 
to highly uncertain primary beam corrections. In order 
of decreasing deep image flux density, these are DR22 
(separation 110'; PB correction -340), DR21 (80'; -25), 
and 4C 44.32 (99'; -108). The last factor listed obeys 
a different Gaussian relation than the flrst two because 
of the different focus settings used in the Kepler and 
Cyg X-3 deep images. 

3.7. Timescales 

We sample the variability of our catalog sources on 
timescales of days to years. If a given source is measured n 
times at dates t^, there are n(n — l)/2 possibly-redundant 
intervals sampled, \ti — tj\ for i < j. We plot these 
intervals in Figure 2 for the three samplings in this work. 
Coverage is fairly uniform across the range of probed 
timescales. 

With additional analysis, our observations can also be 
sensitive to variability on timescales around the observa- 
tion time of a typical epoch image, —1 hour. Evolution at 
and below this threshold can be searched for by imaging 
subsets of the data within each epoch, with the expected 
tradeoff between time resolution and sensitivity. Our 
maximum time resolution is the ATA correlator integra- 
tion time of 10 s. Any single-epoch transients discovered 
in the complete dataset will be examined for intra-epoch 
evolution in this manner, but in this work we flnd no 
sources that merit this detailed investigation. In the gap 
between our —hour integration time and our —week ob- 
serving cadence, we are potentially sensitive to events 
(e.g., a lucky observation may occur precisely during an 
hour-long stellar flare) but can only poorly constrain their 
duration. 

Our sensitivity to slowly- evolving sources might cor- 
respondingly be improved by imaging, e.g., months of 
data at a time, and searching for sources too faint to be 
detected in a single epoch but too variable to be detected 
in the deep image, as in Bower et al. (2007) and other 
works. (This is simply a rough matched-filter approach.) 
Because our images are not limited by thermal noise it 
is unclear how much of a benefit this technique would 
provide in practice. We defer an investigation of this 
approach to future work. 

4. FIRST RESULTS 

In the following subsections we present the first results 
from analyzing the subset of ASGARD data described 
above. 

4.1. Images 

The deep image of the Kepler pointing at 3 GHz is 
shown in Figure 3. All of the observations of this pointing 
were made with a focus setting of 1.9 GHz, so the primary 
beam is broadened as compared to the nominal, optimally- 
focused configuration. 



The optimally-focused deep image of the Cyg X-3 point- 
ing at 3 GHz is shown in Figure 4. The LSS and FOV 
may be compared with the 1.4 GHz image made with the 
Westerbork Synthesis Radio Telescope (WSRT) presented 
in Setia Gunawan et al. (2003, fig. 2; S+03 hereafter). 
The ATA's compact configuration and large FOV make 
it much more sensitive to LSS than the WSRT, and it 
is notable that the image shown in Figure 4 is produced 
from (multiple epochs of) a single pointing. We detect 
significant LSS out to a radius of —45' from the pointing 
center. 

Figure 5 shows an image made from a single Cyg X-3 
epoch (2011 Feb 01) after subtraction of the LSS. The 
Cyg X-3 field contains three bright sources near the half- 
power point that are somewhat problematic in the subtrac- 
tion: DR7 (west), DR15 (south), and 18P61 (northeast; 
Wendker et al. 1991). All of these are at least marginally 
resolved and embedded in cuspy extended emission. The 
residuals due the these sources nonetheless do not signifi- 
cantly impair the imaging of each epoch. Not shown in 
the figures, but easily detectable in the Cyg X-3 epoch 
images, are DR22 and DR21, as described above. 

We generated a total of 29 epoch images, 23 of the 
Cyg X-3 pointing (dropping six analyzed epochs with 
dynamic range issues caused by the 2011 March major 
fiare), and 6 of the Kepler pointing. In Figure 6 we plot 
the representative rms of each image (as reported by sfind) 
as a function of integration time. In Figure 7, we show 
the number of sources detected in each image. 

To assess the effectiveness of the LSS subtraction pro- 
cess, we compare image rms values in regions with varying 
levels of LSS fiux density. We selected several source free 
regions in the Cyg X-3 and Kepler fields of 50x100 image 
pixels. For each epoch image and each region, we com- 
puted the ratio of the rms in that region to the rms in an 
equally-sized region on the outskirts of the image, away 
from all source emission. Taking this ratio compensates 
for the varying noise baseline of each image. Figure 8 
summarizes these ratios as a function of the mean LSS 
model fiux density associated with each source-free region. 
(Measurements from the Kepler field, in which there is 
no LSS subtraction, are assigned a mean LSS model fiux 
density of zero.) The rms ratios increase above the base- 
line as LSS becomes more significant, but saturate at a 
factor of —2 for mean LSS fiux densities of —0.4 mJy. We 
placed a group of regions especially near the three bright 
sources of the Cyg X-3 field. These display somewhat 
higher rms ratios than other regions in areas distant from 
these bright sources, but the difference is slight. 

To give a sense of the areal coverage of the ASGARD 
GC field, we show a preliminary mosaicked image of 3 GHz 
GC data in Figure 9. LSS is not subtracted in this image. 
The missing pointings have coverage at lower frequencies 
but not at 3 GHz. The Sgr A complex, with structured 
emission reaching brightnesses of 50 Jy/beam, presents a 
clear challenge. The data contributing to Figure 9 come 
from only a few epochs, so the more thorough hour angle 
coverage of the complete dataset will make a significant 
difference to image quality. Joint deconvolution of some 
or all of the pointings might greatly improve the deep 
map, although this would depend strongly on how well 
the ATA primary beam can be modeled. 
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4.2. Source Catalog 

We robustly detect 134 compact sources in the deep 
images of the two pointings. There are 86 sources detected 
in the Kepler region and 48 in the Cyg X-3 region, the 
difference between the two being primarily due to the 
effects of LSS in the latter pointing. The faintest source 
discovered in the current Kepler deep image has a flux 
density of 2.6 mJy. The faintest source discovered in the 
Cyg X-3 image has a flux density of 9.2 mJy. 

Because it is 13.5° out of the Galactic plane and con- 
tains no LSS, the Kepler pointing provides a clean testbed 
for our catalog. All but one of the sources associated with 
the Kepler pointing in our catalog has a match to a 
source in the NVSS catalog within 20'^ and none of these 
match multiple NVSS catalog entries. The probability 
of an individual chance match between the catalogs at 
this positional tolerance is ~0.5% (Condon et al. 1998). 
The unmatched source is visible, but marginal, in the 
NVSS imagery. We define a set of NVSS sources that 
might be detected in the Kepler deep image as those 
that are projected to have a primary-beam- attenuated 
(ATA-apparent) flux density greater than 2.6 mJy, so 
long as the PB correction factor /pb is smaller than two. 
Two-thirds (43/66) of these sources are detected, with 
most of the nondetections being marginally discernable 
in the deep ATA imagery. Meanwhile an additional 41 
NVSS sources are detected that do not meet the above 
criteria. The disjunction between the detections is a result 
of some combination of spectral dependence, uncertain 
PB modehng, incompleteness to marginal detections, and 
possible genuine source variability. Of the undetected 
NVSS sources, the maximum predicted ATA flux density 
is 5.6 mJy; that is, there are no NVSS sources that should 
have easily been detected in the ATA data that were not, 
and so there are no bright NVSS sources in the field 
that reduced their flux density by a significant fraction 
by the time of our observations. All NVSS sources with 
predicted ATA flux densities >5.8 mJy and PB correction 
factors less than 2 are detected. If no PB correction limit 
is applied, the limit is 17.7 mJy, with a maximum applied 
PB correction factor of 21.0, corresponding to an angular 
separation of 79^ 

NVSS imaging of the Cyg X-3 region is quite poor and 
so we do not compare the catalogs for this region. Instead 
we use the S+03 L-band catalog of compact sources in 
the Cygnus 0B2 region. Of our catalog entries associated 
with the Cyg X-3 pointing, all but 8 are matched to S+03 
sources. Of those, four are outside of the S+03 FOV, 
one is clearly detected but is too extended to meet their 
selection criteria, and three genuinely do not appear to 
be detected in the S+03 maps, in either their L-band or 
350 MHz observations. Using the same criteria as above, 
one S+03 source might have been expected to be detected 
with a predicted apparent ATA flux density of ~10 mJy, 
but it is blended with the very bright DR15. Ah S+03 
sources with predicted ATA flux densities >11.2 mJy and 
PB correction factors less than 2 are detected. 

In Figure 10 we compare various source counts. In 
order to avoid uncertain primary beam corrections, we 
restrict our assessment to sources within the half-power 
point, reducing our catalog to 43 and 27 sources within 
the Kepler and Cyg X-3 pointings, respectively. With 
such limited numbers, differential functions are extremely 



uncertain, and so we report cumulative source count func- 
tions. In the Kepler pointing, our results are comparable 
to those of the NVSS, which is expected given the results 
reported above. Our Kepler counts are somewhat higher 
than those found in PiGSS-I, probably due to a combi- 
nation of higher-quality deep imaging and manual source 
identification, both made possible by the much smaller 
AS CARD survey footprint. Our source counts in the 
Cyg X-3 pointing are somewhat lower than those of S+03, 
which may reflect the different observing frequencies of 
the two surveys, although such an effect would also be 
relevant to the NVSS comparison. Both the ASGARD 
and S+03 source counts are uncertain on the 20-50% 
level due to small-number statistics. 

4.3. Detection Limits and Completeness 

In order to analyze searches for transient and highly- 
variable sources, we must understand the completeness of 
these searches. The relevant detection limit for transient 
searches such as our sfind step (§3.4) is, by definition, 
the blind source detection limit. Furthermore, because 
the ATA will not resolve true astrophysical transients, 
the detection limit may be expressed as a flux density 
limit rather than a brightness limit. As a shorthand 
we thus refer to this particular quantity as the "BUDL": 
blind unresolved-source detection limit. If artifacts are 
not significant, an image that has not been corrected for 
primary beam attenuation has a single value of the BUDL 
in terms of apparent flux density: the limit reported 
by sfind. Figure 6 depicts the variation of this value 
as a function of integration time for the epoch images 
analyzed in this work. The BUDL in terms of intrinsic 
flux density varies spatially as defined by the primary 
beam response. The apparent BUDL is a function of 
each image's thermal noise; the additional effective noise 
due to limitations in calibration, deconvolution, etc.; and 
the source detection cutoff determined by the sfind FDR 
algorithm. Typically, the apparent BUDL ~ 5.1(7, where 
cr is the representive image rms, but the coefficient varies 
in an image-dependent manner due to the use of the 
FDR algorithm. As shown in Figure 6, however, its range 
of variation is not very large. We generally find that 
BUDL oc T~^l^ as would be expected, where r is the 
image integration time. 

We insert artifical sources into our imagery and then 
apply our source-finding method in order to find an up- 
per limit to our method's completeness as a function of 
flux density. This approach only finds an upper limit 
to the completeness because we insert the sources mid- 
way through our processing pipeline, after the u-v data 
have been initially calibrated. We experimented with 
inserting false sources in both the u-v and image domains 
and obtained equivalent results with both approaches, so 
for more detailed studies we used image-domain inser- 
tion, which avoids the need to rerun the computationally- 
expensive imaging steps. 

We began by sampling the differential completeness 
at fixed apparent flux density CaiSa)-, that is, without 
accounting for primary beam attenuation. To sample 
this function, we spread ~500 unresolved sources of flux 
density Sa throughout all of the epoch images, spacing 
them randomly but evenly in position angle and radius 
relative to the pointing center to avoid overlaps. (This 
is fair because we expect any image to contain at most 
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one transient.) The blind source detection routines are 
then run and the fraction of detected sources is reported 
as the completeness, assigning an uncertainty to the mea- 
surement assuming Poisson statistics. The top curve of 
Figure 11 shows the results of sampling 4 < S'^ < 200 mJy 
and a four-parameter analytic fit to these results. Our 
model function is 

,(c^=i^ iiSa<{-D/BY'C, 
Ca{Oa) |Atanh(B5f + D) otherwise. 

(2) 

There is no particular theoretical motivation for this 
representation, but it can be made to match the data 
well empirically and has two desirable analytical traits. 
Firstly, it plateaus as 5'^ ^ oo. Secondly, it has realistic 
discontinuous behavior at the completeness zero point. 
The fitted parameter A, the plateau value, is 99.3%, and 
reflects that a small fraction of our images is disqualified 
from blind source detection by virtue of proximity to a 
cataloged source. 

We also sample the differential intrinsic completeness 
function Ci{Si), assuming our analytic PB model, by in- 
serting false sources with appropriately PB-attenuated 
flux densities. Our results are also shown in Figure 11, 
along with a corresponding fitted model of the complete- 
ness function in which the plateau value was fixed to 
match that of the apparent completeness function. There 
is substantial disagreement between our samples of this 
function and an estimate derived from Ca and our image 
BUDLs, suggesting that the assumption of a spatially- 
uniform BUDL in each image does not hold strongly. 

Finally, we compute the cumulative intrinsic complete- 
ness Ci] that is, the completeness to all sources intrinsi- 
cally brighter than some limiting value. In order to do this 
we must assume a distribution function for source number 
counts as a function of flux density. We use the standard 
Euclidean, volume-limited distribution N{>S) cx S~^^'^, 
which agrees well with the observed number counts of our 
catalog (Figure 10). The cumulative intrinsic complete- 
ness is then: 

Ci{>Si) = / ds'MSDisl)-'/' / / dslis^r'^' 

JS^ I JSi 

(3) 

We evaluate this integral numerically using the analytic 
fit to our samples of Ci{Si). Below the cutoff of q at 

4.4 mJy, a(>5i) (X^f^ 

4.4. Search for Transients 

We performed a search for transient radio sources with 
our subsample of processed ASGARD images. Recall 
that by our definition, "transient" sources are merely 
ones that are not detected in our deep images, and so 
must be searched for in the individual epoch images. 
Sources detected in this way will, by construction, evolve 
on relatively short timescales, but it is important to 
note that we are also sensitive to sources that vary on 
relatively long timescales: it is just that they are detected 
in the deep images, so they do not satisfy our definition 
of transience. Once detected, any transient sources are 



entered into our catalog and their photometry is measured 
at every epoch, so both classes of sources are analyzed 
identically, although in the stereotypical single-epoch case 
the "lightcurve" for a transient will consist of a numerous 
upper limits and a single detection. 

After applying the sfind source detection and filtering 
procedure described in §3.4, we were left with a list of 17 
transient candidates in our 29 epoch images (23 of the 
Cyg X-3 pointing excluding of the 2011 March major flare, 
six of the Kepler pointing). All but one were detected 
in the Cyg X-3 pointing. Visual inspection confirmed 
that all of the candidates were spurious, as indicated by 
various combinations of unphysical extended structure, 
a large distance from the pointing center, association 
with sidelobes or incompletely subtracted LSS, and/or 
indistinguishability from other noise fluctuations in the 
image. The full-scale ASGARD transient search will verify 
transient candidates rigorously, both by partitioning the 
imaged data to check for instrumental errors (Frail et al. 
2012, cf.) and by treating the image noise statistics 
quantitatively. 

Although experience suggests that systematic effects 
are far more likely to cause false detections than thermal 
fluctuations, it is instructive to consider the limits im- 
posed by noise (Frail et al. 2012). Our search examined a 
total area of ~7 x 10^ synthesized beams, accounting for 
the variable beam size and HPBW of each image. This 
corresponds to one expected noise event of 4.68<j assum- 
ing purely Gaussian statistics. All of our candidates are 
above this threshold, so our search is not yet contami- 
nated by statistical fluctuations. Extrapolating to the 
complete 3 GHz dataset, one statistical noise fluctuation 
of ~5.4<j may be expected. 

We quantify the power of this early search by comput- 
ing its effective search area as a function of the BUDL. 
Given the BUDL of an image Sa lim? our analytic primary 
beam model, and the cumulative apparent-flux- density 
completeness function Ci{Si) determined in the previous 
section, the completeness-corrected effective solid angle 
in which sources with intrinsic flux density greater ^^^um 
may be detected is 

^eff(^a,lim, >5'i,lim) = C'i(>S'i,lim)7rr {Sa,\ im? ^i,lim) 7 

(5) 

where r is the radius within which every source of intrinsic 
flux density >Si will be detected, given a fixed BUDL: 

ro if5z<5a, 
r\Sa, S,) = I 2al^ log(f ) if S, < 905,, (6) 
l(3crpB)^ otherwise. 

(The upper limit stems from our rejection of sfind sources 
found beyond the 98.9% attenuation point.) Following 
Bower et al. (2007), a transient survey consisting of N 
visits to the same field of solid angle O probes a total 
area of r^Q{N — 1), if the time elapsed between visits is 
larger than the transient timescale. If the noise in each 
epoch varies, the effective area searched in each epoch 
also varies, given a fixed S'i,iim- Denoting these areas 11^, 
where > ^2 > • • • > ^at, the effective area searched in 
this case is Xli^2 ^^^^ probed by each epoch is 

a strict subset of the area probed by every more sensitive 
epoch. 

The above condition holds in the case of identical point- 
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ing centers and nonincreasing radial primary beam shapes, 
but not for mosiacs with overlapping pointings. The mo- 
saic case can be treated numerically by expUcitly mapping 
the areal contribution of each epoch. We take this ap- 
proach to determine the effective area as a function of 
'S'i,iim for our processed data. Using the same approach as 
Frail et al. (2012), these measurements can be converted 
into upper limits on snapshot transient source areal den- 
sities. If a total area of r^eff is searched and no transients 
are detected, the 95% confidence limit (CL) density upper 
limit is to a very good approximation. Our results 

are plotted in Figure 12 and discussed in §5. 

We also extrapolate this technique to anticipate the 
parameter space that will be probed by the full survey. 
We use the empirical relationship between integration 
time and BUDL shown in Figure 6 to estimate the limit 
for all epochs as yet unprocessed. When doing so we 
add in a scatter comparable to that seen in the figure. 
This approach "bakes in" typical levels of data loss due 
to RFI, instrumental malfunctions, and so on, assuming 
that the data processed thus far are not unusual in those 
regards. Compared to what is presented in this work, 
the full dataset will be more powerful at the bright/rare 
end because most of the longest integrations have already 
been processed. 

4.5. Variability Analysis 

Several different metrics are commonly used for as- 
sessing the variability of radio sources. Most techniques 
compute a single scalar variability metric for each source, 
although this approach is necessarily reductive. (For in- 
stance, a source may be highly variable on one timescale 
and less so on another.) Analyses such as the structure 
function approach (Simonetti et al. 1985; Emmanoulopou- 
los et al. 2010) are more sophisticated, although they still 
encode certain assumptions about the nature of the vari- 
ability being measured. Well-designed scalar metrics can 
still, however, capture meaningful information about over- 
all variability, as we describe below. 

Becker et al. (2010) define a modulation index or "frac- 
tional variation" 



/ 



5„ 



(7) 



but as discussed by Ofek et al. (2011) this metric has 
irregular statistical properties depending on the number 
of epochs of observations. The same is true of the metric 



V 
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(8) 



used by Gregory & Taylor (1986), which is algebraically 
interchangeable with the above. Ofek et al. (2011) prefer 
the ratio of the standard deviation of the fiux density 
measurements to the weighted mean, written in their 
notation as "StD/(/)" (and sometimes also referred to 
as "the" modulation index). In our notation this would 
be as/{S), where the maximum- likelihood mean fiux 
density (S) of a set of n measurements Si with associated 
uncertainties a,- is 



While as/ (S) is useful for comparing different studies, we 
disfavor it for the ranking of candidate variables because 
it is insensitive to the scale of the uncertainties in each 
measurement, and thus the confidence with which varying 
and nonvarying sources can be distinguished. (This is 
easy to see if the ratio is expressed as 



as ^ V{S^) - (g)- 
(S) 



(S) 



(10) 



where the uncertainties only appear implicitly in the 
averages.) 

Building on PiGSS-II and Bannister et al. (2011a) we 
prefer statistics for this purpose. The statistic 
regarding the hypothesis that some source is unvarying is 



X 



E 



S^ - (S) 



(11) 



(9) 



The distribution of observed values among the ob- 
served sources is not well-defined because, as alluded to 
in Bannister et al. (2011a), they are not drawn from one 
parent distribution. They rather come from the family 
of Xk distributions, where k = n — 1 is the number of 
degrees of freedom associated with each source. Com- 
puting a reduced does not help because that only 
normalizes the expectation values of the distributions, 
not their shapes. To obtain a well-defined distribution we 
must instead apply the full cumulative distribution func- 
tion (CDF) of the Xk family, computing for each source 
the probability of accepting the hypothesis that it is 
constant. We give the expression for Pc in Appendix A 
(Equation Al). In the case of Gaussian errors and no 
variability, the observed Pc values of all sources will be 
uniformly distributed between and 1. 

The metrics listed above do not use the timestamps 
associated with each fiux density measurement. There- 
fore, although they can provide an overall assessment 
of whether a source is somehow "variable," they cannot 
describe the nature of that variability. Among this set of 
metrics Pc is in some sense ideal because it is precisely 
a probabalistic assessment of this matter. Ofek et al. 
(2011) find that the structure function of Galactic radio 
variables saturates at r ~ 10 d, remaining fiat to at least 
r ~ 60 d, suggesting that our basic rankings are useful. 
Nonetheless for particular sources structure functions can, 
for instance, probe the contribution of scintillation to 
the observed variability (e.g., Rickett et al. 2000; Lovell 
et al. 2008). The many epochs of ASGARD provide a 
rich dataset for this analysis: the 3 GHz dataset spans 
a total time baseline of 1.2 yr and contains 18 pointings 
with more than twenty visits and 14 pointings with more 
than fifty. 

Sources farther away from the pointing center tend 
to have increasing systematic effects. This can be seen 
by examining trends in Pc or as / (S) against the PB 
correction factor /pb, as seen in Figure 13: although 
genuine variability will not increase with distance from the 
pointing center, apparent variability does. In Figure 14, 
we compare the Pc CDFs of sources inside and outside of 
the half-power radius. In the latter population, there is an 
overabundance of both apparent variability and sources 
with overestimated uncertainties (cf. Appendix A). This 
leads us to restrict our variability analysis to sources 
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within the half-power point, which constitute 53 of the 97 
sources that have sufficient detections to assess variabihty. 

The excess on the low end of the Pc CDF may be indica- 
tive of the presence of genuinely variable sources, but the 
sample size is insufficient to allow conclusive determina- 
tion. For the same reason we do not attempt to evaluate a 
variability confidence threshold as computed by Bannister 
et al. (2011a). More generally, the irregularity in the CDF 
and its divergence from the theoretically-motivated shape 
defined by Equations A4-A6 indicate that our Pc values 
are not statistically rigorous, although they are useful 
for ranking the variability of sources. For instance, in 
Figure 15 we compare the dependence of Pc and (Js/{S) 
on PB-corrected source flux density. For the more reliable 
sources with /pb < 2, there is a clear increase in as/{S) 
as (S) decreases. As mentioned above, we attribute this 
to the fact that as / (S) does not account for the overall 
scale of the uncertainties in a set of measurements, which 
are fractionally larger for fainter sources. The Pc metric 
does not show this flux density dependence. 

In Table 3 we present positions and variability metrics 
of the eight most variable sources, as ranked by Pc, in 
the subset of sources found within the half-power point. 
Light curves for these sources are plotted in Figure 16, 
and image cutouts showing the sources as imaged are 
shown in Figures 17 and 18. We discuss the nature of 
these sources in §5.2. Cyg X-3 is the most variable source 
and is emphatically detected as so by every metric, even 
without the inclusion of the 20 Jy measurements from its 
2011 March major flare. None of the other sources ranked 
as highly- variable are obviously so. Six of the eight highly- 
variable sources are found in the Kepler pointing, and four 
of these have nearby companions (see §5.2). Although our 
photometry routines simultaneously fit nearby sources, 
it is possible that undeconvolved sidelobes are affecting 
our flux density measurements. Epochs with low-quality 
images tend to be more obviously subject to systematic 
photometric effects. In our investigations, post-imaging 
calibration of these images has been unable to remove 
these effects: although certain sources are systematically 
shifted to lower flux densities, different sources are sys- 
tematically shifted upward, so a simple scale factor is an 
insufficient correction. 

5. DISCUSSION & CONCLUSIONS 

Our first results demonstrate the characteristics of the 
ASGARD dataset, the strategies we have used in our pro- 
cessing pipeline, and an initial search for highly variable 
and transient sources. 

5.1. Galactic Radio Transient Areal Densities 

Figure 12 shows estimates of areal densities for various 
galactic radio transient phenomena and the parameter 
spaces probed by several observational efforts. In this 
section, we describe the different components shown in 
the Figure and their derivations. This standard "logA^- 
log6"' plot inevitably collapses important distinctions 
between different surveys and populations, such as observ- 
ing frequency, relevant timescale, and spatial distribution. 
When relevant we attempt to make these distinctions ex- 
plicit in the discussion below. Most densities on the plot 
are for transient sources, as defined either empirically (for 
the observations) or by an order-of-magnitude increase in 



flux density (for the theoretical predictions). The excep- 
tions are the values shown for extreme scattering events 
(ESEs) and the Becker et al. (2010) variability results, 
both of which correspond to ~50% variability. Although 
we do not estimate their areal density here, we note that 
flares from pre-main-sequence stars are detectable as ra- 
dio transients (Bower et al. 2003; Forbrich et al. 2008; 
Salter et al. 2008) evolving on hour timescales at mm 
and cm wavelengths, and so are an additional potential 
source of Galactic radio dynamism. 

To provide a reference for the values we determine 
below, we also include on Figure 12 the snapshot den- 
sity of extragalactic tidal disruption events similar to 
Swift J1644+57 as estimated by Frail et al. (2012). These 
events are detected at cm wavelengths and evolve over 
timescales of about a month. 

5.1.1. This Work and Forecast for Complete 3 GHz Survey 

We show the parameter space probed in this work and a 
forecast of the results of a complete analysis of the 3 GHz 
survey data as described in §4.4. The forecast comes in the 
form of an areal density upper limit should no transients 
be discovered. This limit should be comparable to the 
extragalactic results of PiGSS-II for sources brighter than 
~10 mJy. Our calibrator pointings could also be used in 
the transient search, but current results suggest the odds 
of a successful detection would be small (Bell et al. 2011; 
Frail et al. 2012). 

5.1.2. M-Dwarf Flares 

We show the snapshot areal density of flaring M dwarfs 
as computed by Osten (2008), 0.11 deg~^, taking her 
"submillijansky" detection limit to be 0.5 mJy. The 
N{>S) (X S~^^'^ scaling here is well-justified because 
M dwarf flares are relatively faint and are only detectable 
out to hundreds of pc at this limit. The rate might be 
increased by similar flares from substellar dwarfs (Berger 
2002). Because coherent flares evolve quickly (r ~ 60 s) 
and can be fairly narrowband (A/// ~ 0.1), the de- 
tectability of these events in practice wih depend strongly 
on specific survey characteristics (Abada-Simon & Aubier 
1997). A survey conducted in the same manner as AS- 
GARD would have extremely poor sensitivity to these 
events, even if it appeared to reach the necessary density 
limit on Figure 12, because of the very different event 
timescales and steep, narrow flare spectra. 

5.1.3. Active Stellar Binaries 

We again follow the analysis of Osten (2008). Drake 
et al. (1989) analyzed VLA observations of 122 RS CVn- 
like active binary systems and detected 66, finding lumi- 
nosity densities ranging from '^10^^-10^^ erg s~^ Hz~^ at 
6 cm wavelengths, with a median of ~10^^ erg s~^ Hz~^. 
We set a cutoff of a factor of >10 luminosity increase for 
such a system to be considered a transient rather than 
a variable. Approximately 10% of the systems observed 
had luminosity densities >10^^ erg s~^ Hz~^. If all of 
these systems are flaring median-luminosity binaries and 
all of the binaries have similar variability patterns, an 
approximate duty cycle for active binaries to appear as 
radio transients is also 10%. We then find that the areal 
density of such systems brighter than 10 mJy should be 
4 X 10~^ deg~^, taking the spatial density of active binary 
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systems to be ~6 x 10~^ pc~^ (Favata et al. 1995). As 
with flaring M dwarfs, these systems are only detectable 
out to hundreds of pc at mJy sensitivities and thus can 
be treated as an isotropic population for our purposes. 
Even with our somewhat generous treatment of both the 
duty cycle and radio luminosity dynamic range of these 
systems, they will be very difficult to detect blindly. Un- 
like M dwarf flares, however, active binary flares evolve 
on day timescales and have flat, broadband spectra that 
are amenable to detection at cm wavelengths (Owen & 
Gibson 1978). 

5.1.4. X-Ray Binaries 

Lutovinov et al. (2007) used the INTEGRAL dataset 
to find 74 XRBs with \b\ < 2°, with 41 of those being 
high-mass systems. Although the distribution of these 
systems is nonuniform and, in the case of high-mass sys- 
tems, appears to be linked to Galactic spiral structure 
(Bodaghee et al. 2012), we compute a characteristic areal 
density in the GP by assuming these sources are dis- 
tributed in a region of |^| < 90°, |6| < 2°. XRBs flare at 
cm wavelengths on hour-to-week timescales, making them 
well-suited to the ASGARD observing strategy (Hjellming 
& Han 1995). We compute typical flaring flux densities 
of these systems as follows. For black hole XRBs, we use 
radio observations of 20 systems presented by Gallo et al. 
(2012). The radio luminosities of systems with multiple 
observations vary by factors of ~10, so that marginally- 
detected flaring systems may appear as transients. The 
maximum observed flux densities of these systems range 
from ~0. 1-400 mJy, with a median of 8 mJy, assum- 
ing flat spectra and an observing frequency of 3.09 GHz. 
For neutron star XRBs, Migliari & Fender (2006) find 
a typical flux density of ~0.4 mJy and a typical radio 
luminosity dynamic range of 5, so that flaring systems 
typically reach ~2 mJy. For both classes of systems, the 
flaring duty cycle appears to be ~1% (Fragos et al. 2009; 
Kording et al. 2005). Combining these results, we ar- 
rive at an areal density of ~10~^ deg~^ for flaring XRBs 
and a characteristic flux density of ~4 mJy, where the 
latter value is intermediate between the ones mentioned 
above with a bias towards a smaller value to reflect the 
smaller luminosity dynamic range of the NSXBs and the 
likelihood that the brighter flaring systems have already 
been discovered. For consistency we plot this estimate 
using S~^^'^ scaling but we warn that this is not as well- 
motivated as in the previous cases. Although XRBs are 
brighter than the aforementioned stellar systems, their 
rarity makes their blind detection difficult. 

5.1.5. Extreme Scattering Events 

Most examples of interstellar scintillation affect source 
flux densities by ~10% (Lovell et al. 2008) and so would 
not be associated with radio transience. Extreme scatter- 
ing events, however, can cause variations of order 50% at 
the ^3 GHz frequencies we consider (Fiedler et al. 1987, 
1994a). These events last months and so are well-suited 
to detection with the ASGARD observing cadence and 
analysis method. Fiedler et al. (1987) find an ESE duty 
cycle of ~7 x 10~^. Considering previous indications of 
increased scintillation in the GP (Becker et al. 2010; Ofek 
et al. 2011) and an association between ESEs and Galac- 
tic structure (Fiedler et al. 1994b; Lazio et al. 2000), we 



assume a doubled ESE duty cycle of 1.4 x 10"^ in the GP, 
which is consistent with pulsar observations (Pen & King 
2012). The areal density of blazars brighter than 100 mJy 
is ~0.03 deg~^ (Padovani et al. 2007), and Kraus et al. 
(2003) detected intraday variability (IDV) in 86% (25/29) 
of a blazar sample they observed. If every source subject 
to IDV may experience an ESE, we find an areal density of 
4 X 10~^ ESE-affected sources intrinsically brighter than 
100 mJy per square degree in the GP. We note again that 
these sources may not be detected as traditionally-defined 
transients because of both the typical scale of the effect 
and the fact that ESEs involve significant dimmings, not 
brightenings, of sources. The S~^^'^ scaling is of course 
appropriate for extragalactic sources, but the rate of ESE 
incidence may vary significantly by line of sight. 

5.1.6. 50% Variables from Becker et al. (2010) 

Becker et al. (2010) defined strong Galactic variables as 
sources with 5'max > 1-56'min (i-e., / > 1.5; Equation 7) 
over the mult i- year time baseline of the observations. This 
is approximately equivalent to as/{S) = 1/3 (Ofek et al. 
2011), where the measurements in question are peak, not 
integrated, flux densities. After correcting for an esti- 
mated extragalactic contribution, Becker et al. (2010) 
find that these sources have an areal density of 1.6 deg~^ 
in their survey, not accounting for the completeness of the 
underlying catalogs. We combine this normalization with 
the CDF of the brightest measurements of each applicable 
source and an approximation of the completeness function 
of the underlying catalogs as determined by White et al. 
(2005), which we assumed applied to both catalogs used 
by Becker et al. (2010). ASGARD observations are not 
directly comparable to those of Becker et al. (2010) for 
two major reasons: firstly, the Kepler pointing that we 
analyze is at a much higher Galactic latitude {b ~ +13.5); 
secondly, our observations lack the 2-15 yr time baseline 
over which much of the variability found by Becker et al. 
(2010) occurred, although the structure function results 
of Ofek & Frail (2011) suggest that the level of variability 
on 1-10 yr timescales is approximately constant. The 
difference between the observing frequencies of the two 
surveys may also be relevant due to the frequency depen- 
dence of scintillation effects in the GHz regime (Hjellming 
& Narayan 1986). Applying the Becker et al. (2010) es- 
timates to our sensitivity limit and the footprint (not 
effective area) of the Cyg X-3 pointing, we would expect 
to detect ~5 highly variable Galactic sources. In our 
catalog Cyg X-3 is the only source with (Js/{S) > 1/3. 
We tentatively attribute this discrepancy to the compara- 
tively short time basehne of our study, but defer a deeper 
analysis until later work. 

5.1.7. GC Radio Transients from Hyman et al (2002, 2006, 

2009) 

In a total of 62 epochs of long-wavelength observations 
of the GC region using the VLA and the GMRT, Hyman 
et al. (2002, 2006, 2009) discovered four robust GC ra- 
dio transients: GCRT J1746-2757, GCRT J1745-3009, 
GCRT J 1742— 3001, and a radio counterpart to the X- 
ray transient XTE J1748-288 (Strohmayer et al. 1998). 
Although this last source was first detected in the X-ray 
band, we count it as an independent radio transient be- 
cause, based on the description of its detection by Hyman 
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et al. (2002), it would have been detected in a blind search 
even without the X-ray detection, and it does not appear 
that the X-ray detection affected the scheduling of the 
radio observations. These observations were performed on 
a ^monthly cadence with several-hour integrations and 
so are sensitive to similar but somewhat longer timescales 
than ASGARD. It is difficult, however, to estimate an 
effective transient search area for this survey because 
the observations were made in heterogeneous conditions 
and the search methodology is not described in detail by 
the authors. Our best estimate of the effective search 
area, attempting to take into account the multiple over- 
lapping pointings, different instrumental configurations, 
and varying sensitivities of each observation, is 120 deg^, 
with a characteristic detection limit of 30 mJy. The full 
ASGARD dataset will probe a significantly larger area at 
this limit. As with the comparison to Becker et al. (2010), 
however, a direct extrapolation of the areal density would 
be inadvisable: Hyman et al. (2002, 2006, 2009) observed 
only toward the GC itself, where the source density is 
larger, and the observing wavelengths of the two searches 
differ by approximately an order of magnitude. 

5.1.8. Galactic Radio Transients from Gregory & Taylor 

(1986) 

Ofek et al. (2011) and Frail et al. (2012) have reported 
a density measurement of 10~^ deg~^ for this survey (in 
their Figures 1 and 6, respectively) but this number is 
erroneous (E. Ofek, 2012, private communication). Gre- 
gory & Taylor (1986) surveyed a footprint of 500 deg^ at 
5 GHz over 16 epochs and discovered one bright (>1 Jy) 
transient, GT 0351+543a, that appeared in one epoch, a 
phenomenology similar to that of the events detected at 
Nasu Observatory (Kuniyoshi et al. 2007; Matsumura et al. 
2007; Niinuma et al. 2007; Kida et al. 2008; Matsumura 
et al. 2009). Both day and year timescales were sampled 
but intermediate ones were not; about 2/3 of the detected 
variables evolved on day timescales. The single transient 
detection straightforwardly gives a density measurement 
of (1.33ti:3o) X 10"^ deg-2 (95% confidence limit), but 
the effective search area of the survey as a function of 
limiting detectable flux density is unclear due to variable 
survey sensitivity. Taylor & Gregory (1983) report that 
over the whole survey the worst-case 3.5cr detection limit, 
corresponding to an empirical 50% survey completeness, 
was 17 mJy. At this flux density limit the effective area 
of the survey is thus no less than 3750 deg^, and a 95% 
CL upper limit on the Galactic transient surface density 
is thus 1.3 X 10"^ deg-^ Taylor & Gregory (1983) also 
report a best-case 3.5cr detection limit of 7 mJy, so that at 
this limit the effective area of the survey is no more than 
3750 deg^, and the Galactic transient surface density at 
that hmit has a 95% CL lower bound of 1.4 x 10"^ deg"^ 
or ~0.6 sky~^. We advise caution in the use of these 
limits because it seems unlikely that the only transient 
detected by Gregory & Taylor (1986) would have a flux 
density >1 Jy, when their survey should have easily been 
sensitive to events <50 times fainter. Gregory & Taylor 
(1986) investigated and deemed unlikely the possibility 
that this event was due to interference, however. No 
source was detected in later VLA followup of the event 
(Tsutsumi et al. 1995). 



5.2. The Nature of the Most-Variable Sources 

We investigated the nature of the most- variable sources 
listed in Table 3 and plotted in Figures 16-18. Here we 
briefly describe our findings, omitting the most variable 
source, Cyg X-3. All of the other sources are consistent 
with being AGN. The sources have values of as/{S) rang- 
ing from 0.10 to 0.18, as compared to 0.69 for Cyg X-3. 
As mentioned in §4.5, none of the sources besides Cyg X-3 
are obviously variable upon visual inspection. 

19:20:15.8 +44:03:05: This source is a marginally- 
resolved double, as conffimed by observations from the 
Cambridge One-Mile telescope (Birkinshaw 1978) and the 
VLA (O'Dea & Owen 1985). Our modehng of this source 
as a single component likely leads to its apparent variabil- 
ity. It is a head-tail radio source found near the center 
of the cluster Abell 2319. The position is consistent with 
that of the infrared source WISE J192015.63+440307.9 
(Cutri et al. 2012), which has colors consistent with an 
AGN (Wright et al. 2010). 

20:30:35.7 +41:06:09: This source is also detected in 
the NVSS (J203032+410634) and the WSRT 327 MHz 
GP survey (Taylor et al. 1996, as 2028+4055). The WSRT 
survey included two epochs of observations separated by 
several years. Taylor et al. (1996) searched for variability 
between the epochs, and this source was not detected as 
a variable. There were no counterparts in other bands 
found in standard catalogs. The flux densities at 0.327, 
1.4, and 3.09 GHz are 150, 12, and 21 mJy, respectively. 
Although variability or an unusual spectrum are possible 
explanations, we interpret the L-band measurement as 
suggesting that much of this source's emission is resolved 
out by the VLA. 

19:23:11.0 +43:50:20: This source is also detected in 
the NVSS (J192310+435016), the Westerbork Northern 
Sky Survey (WENSS, as B1921.6+4344; Rengehnk et al. 
1997), and early Westerbork observations of Abell clusters 
(Harris & Miley 1978). It is 12% hnearly polarized at 
5 GHz (Broten et al. 1986) which may contribute to some 
apparent variability. There were no counterparts in other 
bands found in standard catalogs. 

19:24:28.5 +44:17:08: This source is also detected in 
the NVSS (J192428+441708), WENSS (B1922.9+4411), 
and the Green Bank 6 cm survey (GB6, Gregory et al. 
1996, as J1924+4417). It has a generally decreasing 
spectrum over this frequency range. The position is 
consistent with that of WISE J192428.35+441709.0, which 
has colors consistent with an AGN. 

19:23:06.3 +43:42:38: This source appears to be 
the northwest lobe of a double source with a fainter 
central component. It is detected in the NVSS 
(J 192306+434233). The central component is detected 
in the IR (2MASS J19230940+4341460) and X-rays by 
Einstein (2E 1921.5+4335) and ROSAT (Brinkmann et al. 
1997). These properties suggest that the ASGARD source 
in question is a lobe of an FRII active radio galaxy (Fa- 
naroff & Riley 1974). If that is the case, its flux density 
should be virtually constant. We do not model the faint 
central radio component, which will increase the source's 
apparent variability. 

19:19:35.9 +44:25:03: This source is also detected 
in the NVSS (J191936+442503) and possibly WENSS 
(B1918.0+4419, separation -20''), but not GB6. It has a 
generally decreasing spectrum over the range of frequen- 
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cies in which it is detected. There were no counterparts 
in other bands found in standard catalogs. 

19:22:14.7 +44:11:38: This source is also detected in 
the NVSS (J192214+441137), WENSS (B1920.6+4405), 
GB6 (J1922+4411), and by Harris & Miley (1978). It has 
a strongly decreasing spectrum over this range. There 
were no counterparts in other bands found in standard 
catalogs. 

5.3. Outlook 

There is substantial evidence that the prevalence of ap- 
parent variability increases with proximity to the Galactic 
plane (Spangler et al. 1989; Ghosh & Rau 1992; Gaensler 
& Hunstead 2000; Becker et al. 2010; Ofek & Frail 2011). 
Ofek et al. (2011) suggest that there is also a significant 
increase in the typical variability amplitude, based on 
the difference between their findings from observations 
at \b\ ^ T and those of Becker et al. (2010) at \h\ < 1°. 
Our observations, with \b\ < 0.5°, will allow a strong 
test of this conjecture if the source variability function 
is flat over 1-3 yr timescales. We also probe Galactic 
longitudes from —3° < ^ < 7°, investigating a possible 
increase in radio variability with proximity to the GC 
as has been tentatively measured (Becker et al. 2010) 
and might be expected simply from the increased source 
density towards Sgr A* (Genzel & Townes 1987). 

As demonstrated by Figures 4 and 9, the ASGARD 
dataset allows for sensitive moderate-resolution imaging 
of static large-scale radio structure in the Galactic plane. 
Extended structures in our deep images will include non- 
thermal radio filaments (Yusef-Zadeh et al. 1984; Law et al. 
2008), H II regions (Brogan et al. 2003; Nord et al. 2006), 
and supernova remnants (Gray 1994) at a wide range of 
GC separations with implications for WIMP dark matter 
models motivated by Fermi 7-ray observations (Hooper 
& Goodenough 2011; Linden et al. 2011), the acceleration 
and composition of Galactic cosmic rays (Protheroe et al. 
2008; Crocker et al. 2011), and the energetics of the 
interstellar medium in the GC region (Morris & Serabyn 
1996; Crocker et al. 2010). Although there are formidable 
technical challenges to full polarimetric calibration and 
imaging of the data, these would provide a wealth of 
information on the structure of the GC magnetic field, 
ionization content, and outflows (Nishiyama et al. 2010; 
Law et al. 2011a). 
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APPENDIX 

PROBABILITY OF CONSTANCY FROM A SET OF MEASUREMENTS 

In §4.5 we consider the probability Pc that a source is unvarying, given a set of flux density measurements. In 
this Appendix we give the equation used to compute this value and derive a condition that indicates the presence of 
incompletely controlled systematic measurement errors. 

To compute Pc we must use the full cumulative distribution function (CDF) of the family. The probability of 
accepting the hypothesis that a given source is constant is the probability of finding a value at least as large as the 
one obtained for that source. This is 

Pc = Q{\,^), (Al) 
where Q is the complementary regularized V function 

(A2) 

and in turn T{s^x) is the upper incomplete T function 

V{s,x) = / f-'^e-*dt. (A3) 

J X 

In the theoretical case of no varying sources and truly Gaussian errors, the observed Pc will be uniformly distributed 
between and 1. If these assumptions do not hold — for instance, if systematics are present — the actual distribution 
of Pc values can differ. Denote the CDF of the observed Pc values Fc{p)] that is, Fc{p) is the probability of measuring 
Pc < p for an arbitrary source in the ensemble. We claim that when systematics are controlled, genuine source variability 
can only lead to Fc{p) > p, ignoring variations due to finite sample size. Therefore, data that show Fc{p) < p are 
suggestive of uncontrolled systematics. By construction, 

Fc{0) = 0, 
Fc{l) = 1, and 

dFc/dp > 0. (A4) 
We furthermore argue that in the ideal case the underlying Pc probability density function 

dF 

Up) = ^ (A5) 

must be nonincreasing, because absent systematic measurement errors it is impossible for sources to appear statistically 
less variable (biased towards larger Pc) than ones that are in fact unvarying. Therefore 

d^Fc/dp^ < 0, (A6) 

and combining with Equations A4 we must in fact have Fc{p) > p on the unit interval. As mentioned above, the finite 
number of sources being measured leads to uncertainty in Fc that allows only a probabilistic statement as to whether a 
given observed Fc is consistent with completely controlled systematics. 
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Figure 1. Comparison between flux densities of cataloged sources as measured in the deep images and those measured in individual epochs. 
Agreement is generally very good, especially for those sources for which the primary beam correction is not large. The vertical set of points 
with a large spread corresponds to the highly variable source Cyg X-3. Several sources have a discernable offset between the deep image 
and epoch flux densities; these are seen in Cyg X-3 pointing, and the offset arises from disagreements between the large-radius behavior of 
the optimally focused PB model (used to obtain the deep flux density) and that of the defocused PB model (used to open the epoch flux 
densities). The nominally very bright sources are subject to large, uncertain PB corrections; see §3.6. 
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Figure 3. Deep image of the Kepler pointing, uncorrected for primary beam attenuation. The grayscale is hnear from zero (white) to 
20 mJy/beam (black). The maximum brightness in the image is ~200 mJy/beam. The rms residual after applying the CLEAN algorithm is 
0.5 mJy/beam. The HPBW for this image is 1.27° and is denoted by the dashed circle. The synthesized beam is 65^^ X Al" at a position 
angle of 98°. 
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Figure 4. Deep image of the Cyg X-3 pointing, uncorrected for primary beam attenuation. The grayscale is hnear from zero (white) to 
150 mJy/beam (black). The maximum brightness in the image is ~500 mJy/beam. The rms deconvolution residual is 1 mJy/beam. A 
Gaussian prior image was used in the maximum-entropy deconvolution process, causing the restored image to contain emission on angular 
scales larger than those sampled by the ATA data. The HPBW for this image is 1.1° and is denoted by the dashed circle. The synthesized 
beam is 60'' x 42" at a position angle of 106°. 
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Figure 5. Image of the Cyg X-3 pointing from 2011 Feb 01 with LSS subtracted. Primary beam correction has not been apphed. The 
grayscale is hnear from -5 mJy/beam (white) to 50 mJy/beam (black). (Note that this is different than that used in Figure 4.) The 
maximum brightness in the image is ~380 mJy/beam. The rms deconvolution residual is 0.9 mJy/beam. The synthesized beam is 85^^ x 36^^ 
at a position angle of 83°. 
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Figure 6. Noise statistics for epoch images as a function of on-source integration time t. Triangles: representative rms values as reported 
by sfind. Circles: apparent blind unresolved-source detection limits (BUDLs; §4.3). At fixed t, achieved values will vary with system 
temperature, data flagging, calibration quality, and the source detection cutoff dynamically determined by the sfind FDR algorithm. The 
dashed line shows the characteristic BUDL assuming a scaling of t"-*^/^; for t = 10 min, this corresponds to ~24 mJy. Both sets of numbers 
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Figure 7. Number of sources detected in each epoch image, as a function of the representative image rms reported by sfind. Symbols are 
grouped by pointing. There are 86 cataloged sources associated with the Kepler pointing and 48 with the Cyg X-3 pointing. The detection 
rate for the Kepler pointing is lower because the Kepler deep image is much less limited by systematics than the Cyg X-3 deep image, and 
so contains many more faint sources that are cataloged but cannot easily be detected in the epoch images. 
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Figure 8. Measurements of the increase in image noise as a function of LSS flux density. Source-free regions of 50x100 pixels were selected 
in the Kepler and Cyg X-3 flelds. For each epoch image and each region, the rms was computed and divided by the rms in an equally-sized 
box in the lower left corner of each image, far away from all source emission. Because these values are ratios, they are independent of 
the primary beam correction. Each plotted point gives the mean and standard deviation of these rms ratios over all imaged epochs for a 
particular pointing, source-free region, and focus setting. The abscissa is the mean flux density of the LSS model of that region. Open 
circles denote measurements from regions very near (~20 pixels) any of the three bright sources of the Cyg X-3 field. Areas in which there 
is significant LSS show rms increases of a factor of ~2 versus the source-free baseline. See §4.1 for discussion. 
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Figure 9. Preliminary mosaicked image of ASGARD GC field. Unlike the other maps in this paper, a primary beam correction has 
been applied. The grayscale is linear from -70 mjy/beam (white) to 800 mjy/beam (black). The maximum brightness in the image is 
~3920 mJy/beam. Each pointing has been imaged to a diameter of about 2.2 times the HPBW. Most pointings have only a single epoch of 
observations contributing to the image. 
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Figure 10. Cumulative source counts in this and other surveys. ASGARD measurements are derived from the cumulative distribution 
function of source flux densities as detected in the Kepler and Cyg X-3 deep images, combined with effective search areas based on our 
analytic primary beam models. Reference values are derived from differential source count measurements reported in S+03, the NVSS 
(Condon et al. 1998), and PiGSS— I. The arbitrarily-normalized S~^^'^ line shows the expected scaling for a Euclidean, volume-limited 
distribution. See §4.2 for discussion. 
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Figure 11. Completeness of blind source detection in various regimes, as determined from simulations with the current set of epoch 
images. Apparent completeness measures the detection fraction of sources of a given flux density in images without correction for primary- 
beam attenuation. Intrinsic completeness measures the detection fraction of sources of a single given flux density after correction for PB 
attenuation. Cumulative intrinsic measures the expected detection fraction for all sources brighter than a given intrinsic flux density, 
assuming luminosity function N(>S) oc S~^^'^. See §4.3 for discussion. 
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Figure 12. Measurements of and limits to snapshot areal densities of radio transients. See §5.1 for detailed discussion. All uncertainties 
are to 95% confidence limit. Upper red line: upper limit on the areal density of Galactic transients from this work alone. Lower red line: 
predicted upper limit from the complete 3 GHz ASGARD dataset, should no Galactic transients be found. Blue line: areal density of 
Galactic sources with fractional variability > 50% at 5 GHz on multi-year timescales as measured by Becker et al. (2010). Green lines: 
extragalactic transient areal density limits from PiGSS-II on daily and monthly timescales. See Bower et al. (2011) for more information. 
Points and limit arrows: measurements of Galactic transient rates from the VLA at 5 GHz (Ofek et al. 2011), VLA at 0.33 GHz (Hyman 
et al. 2002, 2006), GMRT at 0.235 GHz (Hyman et al. 2009), and NRAO 91-m transit telescope at 5 GHz (Gregory & Taylor 1986), as 
discussed in the text. Dashed lines: areal density estimates for various radio transient populations as discussed in the text. Horizontal line: 
density corresponding to one event on the whole sky at any given time. 
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Figure 13. Variability metrics as a function of primary beam correction factor /pb- The upper panel plots 1 — Pc so that a higher 
vertical position indicates more variability in both panels. The distributions of the metrics both skew toward higher values as /pb increases, 
suggesting less-reliable measurements. Cyg X-3 is omitted from these plots. It has 1 — Pc ~ 1 and cr/{S) ~ 0.69. 
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^ c 

Figure 14. Cumulative distribution of observed Pc values for sources inside and outside the half-power point: given a randomly-chosen 
source in the dataset, Fc{p) is the probability that its Pc value is less than p. Dash-dotted line: the uniform distribution, Fc{p) = p, which 
would be expected in the absence of variability with purely Gaussian errors. Dotted line: simulated observations based on our sampling and 
a population in which 10% of sources have log- normal flux density variability with a scatter of 0.13 dex. Solid (dashed) line: distribution of 
values inside (outside) the half-power point. See §4.5 and Appendix A for discussion. 
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Figure 15. Variability metrics as a function of deep image flux density. The upper panel plots 1 — Pc so that a higher vertical position 
indicates more variability in both panels. Among the sources with /pb < 2, as/{S) increases for the less reliably-measured faint sources. 
The quantity 1 — Pc does not show an obvious trend. 
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Figure 16. Light curves of the most variable sources in this study. The variabihty metric Pc increases from left to right, top to bottom. 
(Recall that lower values of Pc are associated with higher probabilities of genuine variability.) The most-probably variable source, shown in 
the top-left panel, is Cyg X-3, here identified as 20:32:25.6 +40:57:26. 
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Figure 17. Image cutouts (in standard J2000 RA/dec. coordinates) of the most variable sources in the Cyg X-3 pointing. Time increases 
left to right, top to bottom. The coordinates of each source and the bounds of the white-to-black linear intensity scale are annotated to the 
left of each set of panels. Each panel is 8.2x8.2'. The top panels show Cyg X-3 itself, here identified as 20:32:25.6 +40:57:26. 
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Table 1 

ASGARD pointing centers and 3 GHz summary statistics. 



Field Identifier^ Galactic Coordinates {£, b) # Epochs^ Integ. Time (hr)^ 



gc-07+1 


356.500 


+0.500 


66 


45.29 


gc-07-1 


356.500 


-0.500 





0.00^ 


gc-05+1 


357.500 


+0.500 


66 


46.28 


gc-05-1 


357.500 


-0.500 





0.00^ 


gc-03+1 


358.500 


+0.500 


64 


14.20 


gc-03-1 


358.500 


-0.500 


59 


10.84 


gc-01+1 


359.500 


+0.500 


67 


14.72 


gc-01-1 


359.500 


-0.500 


76 


22.44 


Sgr A* 


359.944 


-0.046 


3 


5.67 


gc+01+1 


0.500 


+0.500 


70 


21.70 


gc+01-1 


0.500 


-0.500 


5 


0.80 


gc+03+1 


1.500 


+0.500 


65 


16.87 


gc+03-1 


1.500 


-0.500 





0.00^ 


gc+05+1 


2.500 


+0.500 


68 


15.92 


gc+05-1 


2.500 


-0.500 





0.00^ 


gc+07+1 


3.500 


+0.500 


33 


4.94 


gc+07-1 


3.500 


-0.500 


67 


8.45 


gc+09+1 


4.500 


+0.500 


70 


19.73 


gc+09-1 


4.500 


-0.500 


26 


2.48 


gc+11+1 


5.500 


+0.500 


69 


11.33 


gc+11-1 


5.500 


-0.500 


3 


6.12 


gc+13+1 


6.500 


+0.500 


72 


16.55 


gc+13-1 


6.500 


-0.500 


18 


14.97 


Kepler 


75.756 


+13.491 


47 


96.72 


Cyg X-3 


79.845 


+0.700 


86 


153.05 


X3+0-2 


79.845 


-0.300 


13 


10.87 


X3+2+0 


80.845 


+0.700 


20 


26.58 


X3+2-2 


80.845 


-0.300 


10 


3.96 



^ Identifiers with numerical codes refer to offsets from a reference point measured in 
half-degree increments. 

^ Summary statistics refer to the 3 GHz subset of data as described in §3.1. 

^ These pointing centers were completely observed by the end of the first GC season, 

before the switch to 3 GHz observing frequencies, and so have no coverage in the 3 GHz 

subset. 



Table 2 

Parameters of ASGARD campaigns. 



Subset^ Field Start Date End Date # Epochs Time On-Source (h) 



Complete 


GC 


2009 May 22 


2009 Oct 14 


64 


186 




GC 


2010 Mar 01 


2010 Oct 21 


84 


285 




GC 


2011 Mar 01 


2011 Apr 11 


22 


54 




Cygnus 


2009 Nov 24 


2011 Apr 11 


159 


377 


3 GHz 


GC 


2010 Apr 29 


2010 Oct 21 


73 


245 




GC 


2011 Mar 01 


2011 Apr 11 


22 


54 




Cygnus 


2010 Feb 03 


2011 Apr 11 


108 


291 


This work 


Cygnus 


2010 Feb 03 


2011 Apr 11 


30 


83 



^ The group of rows labeled "Complete" gives summary statistics for the complete 
ASGARD dataset. The group labeled "3 GHz" gives statistics for the 3 GHz subset of 
data as described in §3.1. The group labeled 'This work" gives statistics for the datasets 
analyzed in this work. 
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+ 43:50:20 
-3 to 32 mJy/bm 



19:24:28.5 
+ 44:17:08 
-3 to 43 mJy/bm 



19:23:06.3 
+ 43:42:38 
-3 to 38 mJy/bm 



19:19:35.9 
+ 44:25:03 
-3 to 18 mJy/bm 



19:22:14.7 
+ 44:11:38 
-3 to 38 mJy/bm 



Figure 18. Image cutouts of the most variable sources in the Kepler pointing. Layout is the same as in Figure 17. 



Table 3 

Eight most variable sources in the ASGARD analysis presented in this 

work. 



R.A. 


Decl. 


(mJy) 


/PB 


Pc 


/ 




20 32 25.6 


+40 57 26^ 


47.5 


1.0 


0.0000^ 


46.39 


0.69 


19 20 15.8 


+44 03 05 


50.2 


1.1 


0.0001 


1.58 


0.13 


20 30 35.7 


+41 06 09 


21.1 


1.4 


0.0036 


1.99 


0.15 


19 23 11.0 


+43 50 20 


25.9 


1.2 


0.0062 


1.54 


0.15 


19 24 28.5 


+44 17 08 


35.8 


1.9 


0.0391 


1.51 


0.13 


19 23 06.3 


+43 42 38 


31.7 


1.4 


0.0432 


1.41 


0.11 


19 19 35.9 


+44 25 03 


15.0 


1.6 


0.0463 


1.69 


0.18 


19 22 14.7 


+44 11 38 


32.0 


1.1 


0.0483 


1.33 


0.10 



Note. — (S) is the weighted mean flux density of the source across 
all epochs, /pb is the primary beam correction factor used to determine 
the intrinsic source flux density. Pc is the probability that the source is 
constant given the measurements, assuming purely Gaussian errors. / is 
the modulation index of Becker et al. (2010). crs/{S) is the modulation 
index of Ofek et al. (2011). See discussion in §5.2. 
^ This source is Cyg X-3. 

^ We compute 7 x lO"^-*^-*^ for this value, but this is certainly subject to 
numerical precision issues. 



